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Abstract A statistical analysis of radial distributions 
of Luminous Red Galaxies (LRGs) from the Sloan Dig- 
ital Sky Survey (SDSS DR7) catalogue within an in- 
terval 0.16 < z < 0.47 is carried out. We found 
that the radial distribution of ~ 106,000 LRGs incor- 
porates a few quasi-periodical components relatively to 
a variable r), dimensionless line-of-sight comoving dis- 
tance calculated for the ACDM cosmological model. 
The most significant peaks of the power spectra are 
obtained for two close periodicities corresponding to 
the spatial comoving scales (135 ± 12) h^ 1 Mpc and 
(101±6) hr 1 Mpc. The latter one is dominant and con- 
sistent with the characteristic scale of the baryon acous- 
tic oscillations. We analyse also the radial distributions 
of two other selected LRG samples: ~ 33,400 bright 
LRGs (-23.2 < M < -21.8) and - 60,300 all LRGs 
within a rectangle region on the sky, and show differ- 
ences of the quasi-periodical features characteristic for 
different samples. Being confirmed the results would 
allow to give preference of the spatial against temporal 
models which could explain the quasi-periodicities dis- 
cussed here. As a caveat we show that estimations of 
the significance levels of the peaks strongly depend on 
a smoothed radial function (trend) as well as charac- 
teristics of random fluctuations. 
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1 Introduction 

It is widely accepted that Luminous Red Galaxies 
(LRGs) are good tracers of the intermediate- and large- 
scale structures of matter in the Universe. The proce- 
dure of spectroscopic identification of LRGs from the 



data of the Sloan Digital Sky Survey (SDSS; e.g. York 
et al.|2000 Abazajian et al.|2009 ) and advantages of its 
employing for statistical investigations were described, 
e.g. by Eisenstein et al. (2001). High intrinsic luminos- 



ity of the LRGs and uniformity of their spectral energy 
distribution allows to identify them at higher redshifts 
than the main galaxy sample (MGS) and so to trace a 
larger volume of the Universe. 

Properties of a spatial large-scale distribution of 
the LRGs were intensively investigated over last years 
(e.g. |Eisenstein et al.||2005[ |Hiitsi||2006| |Percival et al 



2007a|b| |Cabre & Gaztanaga 2009[ |Gaztaiiaga, Cabre 


& Hui 2009; Martinez et aLl[2009l 


Sanchez et al.||2009[ 


Kazin et al. 2010a|b 


Sylos Labini||2010[ |Percival et| 


al. 2010 Blake et al. 


2011) with a special interest 



to a large-scale feature in their two-point correlation 
function and to associated series of features in the 
spherically-averaged power spectrum P(fc) (h~ 3 Mpc 3 ). 
These features are known to interpret as a display of the 
baryon acoustic oscillation (BAO). The majority of the 
cited papers confirm the existence of a significant fea- 
ture in the distribution of matter, although |Sylos Labini| 
et al. ( 2009a[ ) came to an opposite conclusion from the 
analysis of the SDSS Main Galaxy Sample (MGS) - 



Kazin et al. 



data release 7 (DR7; see also |Sylos Labini, Vasilyev fc| 
Baryshev||2009b I 



2010a I calculated the two-point corre- 
lation function using a sample of LRGs from the SDSS 
DR7. In the framework of the ACDM cosmological 
model the authors have obtained the baryon acous- 
tic peak at (1—2) a significance level and determined 
the peak position as (101.7 ± 3.0) hr 1 Mpc within the 
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Fig. 1 Angular distribution of LRGs from SDSS DR7 over 
the sky in the Equatorial coordinates; all grey coloured regions 
comprise DRl-Full sample of data (105 831 LRGs) represented 
in Table 1 by |Kazin et al.| l |2010a[ |; the data are available on the 
World Wide Web (see text). 



cosmological redshift interval z = 0.16 — 0.36. They 
used the mock galaxy catalogues produced by the Large 
Suite of Dark Matter Simulations (LasDamas^jto esti- 
mate a sample variance and systematic errors of the cal- 
culations and thus to make their results more reliable. 



This result has been confirmed recently by Blake et al. 
(2011) for an extended redshift range z = 0.16 — 0.44 
with a greater significance (3.4 a) and a slightly spec- 
ified position, (102.2 ± 2.8) h' 1 Mpc, of the baryon 
acoustic peak. 

On the other hand 



Ryabinkov & Kaminker (2011) 



have shown that the radial (line-of-sight) comoving 
distribution of absorption-line systems (ALSs) regis- 
tered in quasar (QSO) spectra within the interval z — 
0.0 - 4.3 reveals a period of (108 ± 6) hr 1 Mpc or (al- 
ternatively) a temporal interval (350 ±20) h~ x Myr for 
the standard cosmological model. The proximity of the 
scales obtained for the spatial LRG distribution, e.g. by 
and the distribution of ALSs by 



Kazin et al. 



(2010a 



Ryabinkov & Kaminker (2011) suggests to analyse the 



radial distribution of the same LRG sample employing 
the technics used by the latter authors. 

In this paper we deal with the SDSS DR7 LRG sam- 



ple b y|Kazin et al.| ( |2010a[ ) presented on the World Wide 
Web|^J The interval of cosmological redshifts consid- 
ered here, z = 0.16 - 0.47, corresponds to the DR7-Full 
sample characterized in their Table 1. The SDSS LRG 
regions on the sky are shown in Fig. [I] 

The basic quantity of the present study is a ra- 
dial distribution function N(rj) integrated over angles a 



(right ascension) and 6 (declination), which is an ana- 
logue of the comoving number density n(z) in the red- 
shift space (e.g. |Zehaviet aL|20 05); T)(z) is a dimension- 
less line-of-site (radial) comoving distance between the 
observer and galaxies under study, N(ji)dr) is a number 
of LRGs inside an interval dr/. 

The dimcnsionless radial comoving distances are cal- 



culated according to the equation (e.g. 


Harrison 


1993 


Kayser, Helbig & Schramm||1997[ Hogg 


1999 


1: 



\Mn(l + 20 3 +^A 



dz, 



(1) 



where i is a numeration of all LRGs, i = 1,2, ...,N to tj 
e.g., N tot = 105831 for the DR7-Full sample by |Kazin 



et al. (2010a I. Following Kazin et al. (2010a I we use 



the dimensionless density parameters O m = 0.25 and 
tt\ = 1 — f2 m = 0.75. The corresponding line-of-sight 
comoving distances are D c .i = c- rji/Ho = 2998 • rji ■ 
Mpc, where Ho = 100 h km s _1 Mpc" 1 is the 
present Hubble constant, c is the speed of light. 

We focus our attention on the search of periodicities 
incorporated in the radial LRGs distribution in rela- 
tion to the ^-variable. We shall treat them as quasi- 
periodicities meaning limited intervals of line-of-sight 
distances under consideration, as well as variations of a 
peak positions and amplitudes for different samples of 
LRGs. In Section 2 we examine power spectra calcu- 
lated for two samples of LRGs employing a point-like 
statistical approach: the full sample DR7-Full and a 
sample from the rectangle region described below. The 
point-like approach is usually used in order to avoid 
sensitivity of results to the procedure of a smoothed 
distribution function (trend) subtraction and to the re- 
lated choice of an averaging bin. For comparison in 
Section 3 we discuss power spectra calculated in a bin- 
ning approach for the DR7-Bright sample which is also 
represented in Table 1 by Kazin et al. (2010a). In Sec- 



tion 4 we use two mock galaxy LasDamas catalogues 
for more conservative estimations of significance lev- 
els of the power spectra, and outline shortly another 
(widely discussed) way to estimate the significance of 
the maximal peaks in the power spectra. Conclusions 
and discussions of the results are given in Section 5. 
Ambiguity effects of a trend elimination procedure on 
results of Fourier spectra analysis is shortly outlined in 
Appendix. 



2 Quasi-periodicity of the ^-distribution: 
point-like approach 



1 http: / /lss. phy.vanderbilt.edu/lasdamas/mocks. html 
2 http://cosmo.nyu.edu/~eak306/SDSS-LRG.html 



Fig. [2] demonstrates the radial distribution function 
N(?7) calculated for independent bins with a narrow 
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Fig. 2 Radial distribution function N(rj) of 105 831 LRGs 
belonging to the grey regions in Fig.[l]and a redshift interval 
0.16 < z < 0.47 or dimensionless radial distances 0.16 < r\ < 
0.43; N(r?) is calculated with using independent bins 8 V — 
0.00033 or c • Srj/Ho = 1 h~ x Mpc; thick solid line displays 
a trend Nt r (f?) determined as a set of parabolas (see text). 



width S v = 0.00033 or 



c • Sn/Ho = 1 hr l Mpc at the 
scale of comoving distances. Narrow spike-like varia- 
tions of N(?7) are noticeable on the background of the 
trend N tr (?7) drawn by a thick solid line. The function 
N tr (r?) can be treated as a smoothed function filtering 



out the largest scales (e.g. Zehavi et al. 2005 Kazin 



et al. 2010a). We assume that the trend can be de- 



termined by selection effects and/or large-scale fluctu- 



ations in the distribution of galaxies (e.g. Sylos Labini 
. In Fig. [2] the function N tr (^) is calculated by 
the least-squares method using a set of parabolas as a 
regression function for N(?/). 

Our aim is to reveal weak but quite significant quasi- 
periodical components and separate them from the rest 
determinate part of the radial distribution (trend), as 
well as fluctuations. However, in the case of such a com- 
plex trend as displayed in Fig. [2] the procedure of the 
trend subtraction is ambiguous in principle and may 
bring to wide variations of peak positions and ampli- 
tudes in power spectra, as demonstrated in Appendix. 
So it is more safely in that case to employ an out-of-bin 
(point-like) statistical technique, which is not sensitive 
to the procedure of trend determination. 

To verify the periodicity of the ^-distribution we cal- 
culate a power spectrum for the sequence of LRG points 
rj{zi) using the Rayleigh power: (e.g. Brazier 



V: 



1994): 



P(m) = 



1 

"Ntot 



/ 27rm?7j 



'N tot 

E cos v 1 



E. / 27rmr/ i 



(2) 



m is an integer harmonic number, L, ; = rj max — r\ m in 
is an interval under consideration. A periodicity yields 
a peak in the power spectrum "P=P(m), with the con- 
fidence probability or cumulative distribution function 
(e.g. |Scargle|p82l ) 



P = [1 - exp(-P)], 



(3) 



where j3 is defined relative to the hypothesis of the Pois- 
son distribution of rj,. This estimation is valid for a 
single independent peak at arbitrary m and yields the 
probability of pure noise generating a power P(m) less 
than given level V . 

Fig. [^represents power-spectra P(fc) calculated with 
the use of Eq. j2j at 1 < m < 40, where the wave 
number k (h Mpc" 1 ) substitutes for m according to 



an expression k — 27rm/D^, where 



is the whole comoving interval. The upper and lower 
panels correspond to two overlapping samples of LRGs, 
both of them belong to the same interval 0.16 < 77 < 
0.43 (L, = 0.27, = 809.4 /i" 1 Mpc). The upper 
panel is calculated for the whole region marked in Fig. [I] 
comprising the main sample of 105,831 LRGs, the lower 
one - for a rectangle region within the central oval-like 
domain restricted by the intervals of right ascension 
140° < a < 230° and declination < 6 < 60°. The 



region is chosen by analogy with Sylos Labini (20101 to 
minimize possible effects of the irregular edges of the 
central domain in Fig. [I] The latter sample contains 
60,308 LRGs. 

The most significant peaks of the spectra P(fc) are 
noticeable for two close periodicities at k = 0.047 and 
0.062 h Mpc -1 or harmonic numbers m=6 and m=8. 
These values k correspond to dimensionless radial co- 
moving scales of rj A,, = 0.045±0.004 and 0.034±0.002, 
or the comoving spatial scales AD C = (135 ± 12) and 
(101±6) h~ x Mpc, respectively. The latter peak (m=8) 
is dominant for both samples with amplitudes well ex- 
ceeding the level 5a estimated with Eq. ([3]). While the 
first one (m=6) is somewhat less significant being lower 
than the level 5a for the rectangle region. 

Fig. [4] demonstrates a reciprocal Fourier transform 
reconstructed for 5 harmonics from m=5 till m=9 and 
shown as a function of the comoving radial distance D c . 
The appropriate direct Fourier transform including am- 
plitudes and phases of harmonics is calculated for the 
sample of all LRGs from the grey regions in Fig. [l] 
One can see that the pair of close periodical compo- 
nents m=6 and m=8 yields a resultant quasi-periodical 
fluctuations of the radial distribution of LRGs. An av- 
eraged scale of the quasi-periodical fluctuations turns 
out to be (100 ± 8) h' 1 Mpc, i.e. it is rather close to 
the main tone m=8. 
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Fig. 3 Power spectra P(k) calculated in the point-like 
approach according to Eq. ([2| for the redshift interval 
0.16 < z < 0.47 (0.16 < T) < 0.43), it is a wave number; 
the horizontal dash lines specify the significance levels 3a, 
4a, and 5<r estimated with the use of Eq. |3| at fixed levels 
of the confidence probability f3 = 0.998, 0.99994, 0.9999994, 
respectively. Upper panel: P(k) calculated for all LRGs from 
the grey regions in Fig. [I] Lower panel: P(fc) calculated for 
LRGs from a rectangle region inside the central oval-like 
domain in Fig. [I] restricted by intervals of right ascension, 
140° < a < 230°, and decimation, < 8 < 60°. 

Less significant peaks are also visible in Fig. [3j A 
wide peak (significance 4tr) at k = 0.217 h Mpc -1 
(m=28), which is present in the both spectra, and a 
narrow peak at k = 0.124 h Mpc -1 (m=16) appearing 
only in the lower panel. The first period corresponds to 
the spatial scale AD C = (29 j^) Mpc. The second 
one, (51 ± 2) hT x Mpc, turns out to be very close to 
the second harmonic of the main peak (m=8). Note a 
difference of the whole sets of spectral features for both 
samples under discussion. 

Although the point-like approach avoids the proce- 
dure of a trend determination a steep rise of both the 
curves P(fc) at the lowest k < 0.05 h Mpc -1 may be 
attributed to the effects of large-scale fluctuations in- 
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Fig. 4 Reciprocal Fourier transformation for five har- 
monic numbers m (see Eq. from the interval 5 < m < 9 
as a function of the comoving line-of-sight distance D c . Ver- 
tical dash lines indicate positions of maxima of the recon- 
structed distribution. An averaged scale between neighbour 
dash lines is (100 ± 8) h' 1 Mpc. 

corporated in the smoothed 77-dependence (see Fig. [2]). 
On the other hand, the long-wave parts of the spectra 
are different for the two samples considered in Fig. [3j 
Such a difference appears, in particular, in relative de- 
crease of the peak at k = 0.047 h Mpc -1 and can- 
cellation of spectral components within an interval of 
0.03 < k < 0.04 h Mpc- 1 in the lower panel. The 
difference could be associated with effects of irregular 
edges of the oval- like region and a set of strips in Fig. [I] 
excluded in the rectangle region. 

Our special simulations show that the point-like 
(out-of-bin) approach is more sensitive to the presence 
of periodicities in a distribution of points than the bin- 
ning approach with relatively wide bins. The latter 
one needs an appropriate procedure of a trend subtrac- 
tion. In the case of simple trends, e.g. close to linear 
or parabolic dependencies, amplitudes and positions of 
peaks in the power spectra calculated in the binning 
mode converge to the values found in the point-like 
approach at successive reducing of a bin width. This 
convergence can be demonstrated, e.g. for the narrow 
independent bins 6 n — 0.00033 used in Fig. [2j How- 
ever, we found that in many cases both the convergent 
procedures yield overestimated significance of peaks in 
the power spectra. It suggests to re-estimate the signif- 
icance levels of the spectral peaks in Fig. [3] in a more 
robust way than it can be produced with using Eq. (§. 
We realize this suggestion in Section [4] 
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3 Bright LRGs: binning approach 

In the cases of quite simple trends (linear or parabolic 
functions) it is more relevant to use the binning ap- 
proach and a proper procedure of a trend, N tr (?7), sub- 
traction. For instance, one can calculate so-called nor- 
malized radial distribution function: 



NN(? ? ) = 



Nfa) - Ntrfo) 



(4) 



Note that instead of the radial distribution func- 
tion N(rj) one can use a comoving number density 
71(77) = N(jj)/(dV/d»7), where dV/d77 is a comoving 
differential volume, which is a simple variation of the 
more conventional value n(z) (e.g., |Zehavi et al.||2005 



Kazin et al. 2010a). This replacement does not change 
Eq. (|4|) and results of following calculations because 
in that case one should divide into the comoving vol- 
ume both numerator and denominator (remind that 
<j(n) = <j(iV)/V, where a is the mean squared devi- 
ation) . 

An example of quite simple trend is demonstrated in 
the upper panel of Fig. [5j We explore the DR7-Bright 



sample indicated in Table 1 by Kazin et al. (2010a 



which contains 33,356 SDSS "Bright" LRGs (-23.2 < 
M < —21.8, M is the absolute magnitude) within a red- 
shift interval 0.16 < z < 0.44 (0.155 < 77 < 0.40). The 
upper panel in Fig. [5] displays the function N(r)) calcu- 
lated for the same independent bins (5 n — 0.00033) as 
in Fig. [2] The trend N tr (77) is calculated as a parabolic 
function by the least-squares method and drawn by a 
thick solid line. 

The middle panel in Fig.[5]represents the power spec- 
trum obtained for NN(7y) according to the equation: 



Pb(fc) 



1 



gNN.sinf 2 ^) 

3=1 V " J 



f 27rm77j 



(5) 



where At> is a number of bins along 77-axis, r\j is a loca- 
tion of bin centers, r]j±d ri /2; j = 1, 2, At, numerates 
bins, in the case of the middle panel At> = 736, — 
Vmax ~ Vmin = 0.245 (or D|; = 735 /i -1 Mpc) is the 
whole interval. The power spectrum contains two sig- 
nificant peaks exceeding the level 3cr at k — 0.017 and 
0.120 h Mpc -1 . The first peak (m=2) corresponds to 
a half of the whole interval and may be interpreted as 
a residual large-scale fluctuation which could be incor- 
porated into the trend Nt r (7?) but it is missed out by 
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Fig. 5 Upper panel: Radial distribution function N(?7) 
calculated for a sample of SDSS "Bright" LRGs (-23.2 < 
M < -21.8) within a redshift interval 0.16 < z < 0.44 
(0.155 < 77 < 0.4) indicated as DR7-Bright in Table 1 by 
|Kazin et al . (2010a). N(t;) is calculated with using the same 
independent bins (S v = 0.00033) as in Fig. [2] solid thick 
line displays a smoothed function (trend) Nt r (7j) drawn as 
2 d -degree polynomial. The whole sample contains 33,356 
LRGs. The vertical dash line shows a shifted low bound- 
ary rj for a reduced sample of the "Bright" LRGs. Middle 
panel: Power spectrum Pb(fc) of the function NN(7?) cal- 
culated according to Eq. (|5| for the whole sample. Lower 
panel: Similar power spectrum Pb(fc) but calculated for the 
reduced sample of 32,395 LRGs with shifted low boundary: 
0.172 < 77 < 0.4. The horizontal dash lines indicate the 
same significance levels as in Fig. ^ 

the parabolic approximation. The second peak corre- 
sponds to a period of AD C = (52.5 ± 2) h^ 1 Mpc, and 
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may be compared with the peak k = 0.124 h Mpc -1 in 
the lower panel of Fig. [3j 

The lower panel in Fig. [5] shows a complementary 
power spectrum calculated in the same way but for a 
somewhat diminished sample of 32,395 "Bright" LRGs 
within a reduced redshift interval 0.18 < z < 0.44 
(0.172 < 7/ < 0.40), i.e. lower edge is shifted. Thus in 
calculations with Eq. |5| we use the values Afb = 685, 

= fjmax ~ Vmin = 0.228. Similar to the middle panel 
the power spectrum contains two significant peaks at 
k = 0.018 and 0.120 h Mpc -1 . The first peak, slightly 
shifted in respect of the peak in the middle panel, im- 
plies the same interpretation. The second one displays 
the same period as the full sample of "Bright" LRGs 
but with higher amplitude ( > 4a). The difference 
between two spectra may be treated as an edge effect 
which provides better tuning of the Fourier component 
at k = 0.12 h Mpc -1 in the case of reduced sample. 

Comparing Figs. [5] and [3] one can notice that the 
power spectra are strongly dependent on a sample of 
objects under investigation. In particular, the radial 
distribution of the bright LRGs displays the peak struc- 
ture of the power spectrum which differs from that ob- 
tained for the whole sample. 



4 Significance estimations: LasDamas 
catalogues 

In this Section we employ two mock galaxy LasDamas 
(LD) catalogs "lrgFull-real" and "Irg21p8-real'j^] (the 
latter simulates SDSS sample of bright LRGs with 
M < —21.8) for complementary estimations of the sig- 
nificance of spectral peaks appearing in Figs. [3] and [5j 
In both cases we use the data of 80 "ns" (North-South) 
realizations. Our aim is to take into account possible 
clustering of the galaxy distribution incorporated in the 



mock LD catalogs (e.g. Berlind & Weinberg 2002 see 



also Manera et al. 2012 ) in contrast to the hypothe- 



sis of Gaussian fluctuations used in Sections [2] and [H 
We calculate power spectra for all realizations of both 
LD catalogs and approximate the whole array of spec- 
tra by the unified ^-distribution applied to spectral 
amplitudes at each harmonic number m (wave number 
k) under investigation. Quality of the a approximation 
for each m or k is controlled by the Kolmogorov criteria 
with a tabulated value Ai_ p at fixed significance level 
p = 0.2 (Ao.g = 1-07). 

To obtain two sets of the power spectra we perform 
calculations of the radial distributions of the LD galax- 
ies with using of the binning approach. We employ the 

3 http: / /lss. phy.vanderbilt.edu/lasdamas/mocks/gamma 
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Fig. 6 Power spectra for two samples of LRGs discussed in 
Sections[2]and[3]in comparison with appropriate significance 
levels (dash lines) and mean power spectra (thin solid lines) 
calculated on the basis of LasDamas mock catalogs (see text 
for details). Upper panel: same power spectrum P(k) as in 
the upper panel of Fig. [3] significance levels of 2a, 3a and 
4cr, as well as the mean spectral line, are calculated with the 
use of 80 realizations of LasDamas "lrgFull-real" catalog. 
Lower panel: same power spectrum Pb(fc) as in the middle 
panel of Fig. [5] significance levels of 2a and 3a, as well as 
the mean spectral line, are calculated with the use of 80 
realizations of LasDamas "Irg21p8-real" catalog. 



same bin S v = 0.00033 as in Section [3j As it was dis- 
cussed in Section [2] the power spectra calculated with 
such a small bin in the case of a simple trend display in- 
significant difference with the power spectra obtained 
in the point-like approach. This allows one to apply 
significance levels estimated on the base of LD data to 
the point-like statistics used in Section [2j 

To make statistical properties of the mock samples 
of the "lrgFull-real" catalog more comparable with the 
main sample of LRGs used in Section [2] we implement 
a reduction of the data, which reconciles simple radial 
selection functions (trends) N^ r D (?7) of the LD data with 
the complex trend of the LRG sample. The procedure 
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of the reduction is performed for all realizations of the 
LD catalog within the interval 0.16 < r\ < 0.40 accord- 
ing to the formula: 



t LD/ 



(0) 



where r\j is a centre of j-th independent bin, Ni n (r]j) 
and Nf in (r]j) are initial and final radial distributions 
of mock galaxies over all independent bins, N^. DSS (?7 J ) 
is a trend of the radial distribution calculated for the 
main SDSS sample, N|~(?7j) is a linear trend calcu- 
lated for each mock realization. For the whole sample 
of 105,831 LRGs we calculate a trend N^ DSS (ry :( ) as a 
reciprocal Fourier transformation of the first five har- 
monics (m=l, 5) with their amplitudes and phases 
obtained in the direct Fourier transform. This trend is 
shown in Fig. [8] of Appendix. 

Then the normalized radial distributions NNj can 
be determined by Eq. Q, where ^finiVj) an d a mean 
value Nfi n stand for N(?7) and N tr (f?), respectively. The 
mean N/ in over the interval 0.16 < rj < 0.40 is used 
as the simplest trend to bring the binning approach of 
Eq. ([5] ) closer to the point-like approach of Eq. pj. 
It allows us to compare statistically a resultant sample 
of the mock power spectra obtained for all realizations 
of the "lrgFull-real" catalog with the spectrum in the 
upper panel of Fig. [3j 

For the samples of the "Irg21p8-real" mock catalog 
we also employ the procedure described in Section [3] 
with the use of Nj„(?7j) and the parabolic trend N tr (i7j) 
for calculations of the values NNj. Following Eq. ([5| 
we obtain a set of power spectra and compare the dis- 
tributions of their amplitudes with the power spectrum 
in the middle panel of Fig. [5| 

It was found with the use of the Kolmogorov crite- 
rion that the ^-distribution f peak amplitudes P(fc) 
is valid for a majority (but not all) of k in the power 
spectra obtained for "lrgFull-real" and "Irg21p8-real" 
realizations. Specifying the confidence probability /3 
we calculate the significance levels of P(k) at each 
k = 27rm/D^, determined by all integer m m i„ < m < 
40; it is performed within the interval 0.04 < k < 
0.25 h Mpc" 1 (m mm = 5) for the whole SDSS LRG 
sample and at 0.017 < k < 0.25 h Mpc" 1 (m min = 2) 
for the "Bright" LRGs. 

Results of the statistical estimations of the peaks sig- 
nificance for both SDSS LRG samples are represented 
in Fig. [6j The upper panel displays the same power 
spectrum as in the upper panel of Fig. [3] and the lower 
panel corresponds to the middle panel of Fig. [5j Thin 
solid lines in both panels are drawn through arithmeti- 
cal means of peak amplitudes calculated for each m 
(m mi „ < m < 40) or respective k for the appropri- 
ate sets of power spectra. Dashed lines in both panels 



are drawn through significance levels 2cr, 3cr, and 4<r 
(upper panel) calculated as the amplitudes P(k) at all 
respective k matching with the confidence probabili- 
ties of the x 2 -distributions: /3 = 0.95, (3 = 0.998, and 
/3 = 0.99994, respectively. Let us note that our model 
estimations are quite approximate due to lack of statis- 
tics of the mock catalogs and may underestimate real 
significance of the peaks. 

In the upper panel of Fig. [6] one can see only two 
peaks k = 0.047 and 0.062 h Mpc -1 exceeding the sig- 
nificance level 3cr, the latter one is higher than 4cr. The 
significance of the third peak k = 0.217 h Mpc^ 1 hardly 
approaches to the level 3er in contrast with its level (4er) 
in Fig. [3j Three noticeable peaks in the lower panel of 
Fig. [6] turn out to be less significant ( > 2a). However, 
the growth of the peak at k = 0.12 h Mpc -1 in the 
lower panel of Fig. [5] as well as an appearance of the 
peak at k rj 0.12 h Mpc -1 in the lower panel of Fig. [3] 
give an evidence in favour of possible presence of the 
additional comoving scale AD C = (52.5 ± 4) hr 1 Mpc 
in the radial distribution of LRGs. This peak should 
be confirmed or rejected by further investigations with 
extended statistics. 

It was marked in Section 2 that the estimations of 
the confidence probability according to Eq. ^ imply 
appearance of a single peak at some wave number k 
in a power spectrum with an amplitude V . Following 



Scargle ( 1982 ) (see also |Frescura, Engelbrecht, Frank 



2008 and references therein) one can consider a set 
of many independent wave numbers fc m ; m= 1, N^, 
where may be ^ 40 used in Fig. [3j and treat any 
of power peaks P(m) as a result of Gaussian noise. 
Then one can estimate so called false alarm probabil- 
ity, F = Pi-(P ma x >V) = l- P Nk , where /3 is defined 
in Eq. (J3|, i.e. probability of at least one of peaks P max 
being equal to (or above) a maximal level V . 

Considering a set of natural wave numbers k cxm for 
the whole sample of LRGs in Figs. [3] and [6] we regard 
that TV*, is close to the theoretical value Nb /2 ( Frescura, 



Engelbrecht, Frank 2008) and exploit also the proxim- 



ity of the point-like and binning approaches at the same 
S 7] = 0.00033. One can deal with the normalized power 
amplitude P„ laT = V = 26.4 at k= 0.062 and u se the 
formula V = -ln[l - (1 -p ) 1 /Nk ] dScargle|l982 



Fres- 



cura, Engelbrecht, Frank|2008 ) ) to match a given level 



of the false alarm probability po with an amplitude Vq. 
In our case Nh = 818, Nk = 409 and we obtain that the 
confidence level (1 — po) = 0.9999994 (significance 5a) 
corresponds to Vq — 20.3. Thus the obtained value lies 
between the 5cr-level in Fig. [3] ( < 15) and the ampli- 
tude 26.4 mentioned above. The estimation makes the 
significance of the main peak lower but not so much 
lower as the 4<r-level in Fig. [6] 
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Analogous estimations can be carried out with the 
power peak at k = 0.12 in the middle panel of Fig. [5j 
The confidence probability of the peak 1 — po = 0.986 
also becomes lower (significance < 3er). In general, the 
estimations employing the false alarm probability yield 
systematically lower confidence levels of the maximal 
power peaks, but these levels are higher than the es- 
timations on the base of LD mock catalogs discussed 
above. Therefore we treat the single-peak estimations 
with using Eq. ^ and the LD-estimations as maximal 
and minimal boundaries of possible confidence levels. 



5 Conclusions and discussion 

The main conclusions of the statistical analysis of radial 
(line-of-sight) distributions of the LRGs within the red- 
shift interval 0.16 < z < 0.47 or the interval of dimen- 
sionless line-of-sight comoving distance 0.16 < rj < 0.43 
can be summarized as follows: 

(1) The radial distribution of - 106,000 LRGs in- 
corporates a few (< 4) significant quasi-periodical com- 
ponents in relation to the smoothed function (trend). 
The most significant peaks of the power spectra are 
displayed for two close periodicities corresponding to 
the spatial characteristic scales (135 ± 12) and (101 ± 
6) hr 1 Mpc. The latter scale is a dominant with signif- 
icance > 4cr. This pair of close periodical components 
may be treated as an integrated quasi-periodical scale 
(100 ±8) h' 1 Mpc. 

(2) Once more appreciable quasi-periodical scale, 
(29 Ij) Mpc, arises in the power spectra as a wide 
peak at a level > 3<r. The peak displays similar am- 
plitudes for two overlapping samples: the main sample 
of 105,831 LRGs and a sample of 60,308 LRGs from 
the rectangle region (see Sect. [2J). Still less reliable 
scale (52.5 ± 4.0) hr 1 Mpc appears for the sample of 
LRGs from the rectangle-region and for the sample of 
"Bright" LRGs. These periodicities must be verified by 
further examinations. 

(3) The main quasi-periodical scale found in this 
work ((101 ±6) h~ Y Mpc) is consistent with the scale 
(102.2 ± 2.8) hr 1 Mpc of the baryon acoustic peak in 
the two-point (monopole) correlation function of the 
LRG distribution (Bl ake et al.||2011 ). These scales are 
in agreement also with a period of (108 ±6) h^ 1 Mpc 
revealed by Ryabinkov & Kaminker ( |2011 ) in the radial 
distribution of QSO ALSs for a wider redshift interval 
0.0 < z < 4.3. Note that the nearness of the periods 
~ lOO/i -1 Mpc revealed at redshifts z < 0.5 for LRGs 
and at z < 4.3 for QSO ALSs argues for spreading of 
the same large-scale periodicity over long interval of the 
cosmological evolution (see also Demiaiiski et al.|2011 ) . 



(4) A set of peaks in the power spectra of Figs. [3j [5] 
and their variations from sample to sample probably ev- 
idence in favour of rather spatial than temporal nature 
of the quasi-periodicities discussed here. Our model 
simulations performed for partly ordered structures of 
points (e.g., Ryabinkov & Kaminker 2011) show that 
it is possible to get a point-like power spectra similar 
to those represented in Figs. [3] and [5] performing simu- 
lations of a cloud-like 3D-distribution of points around 
vertices, e.g. of a face-centered cubic. 

(5) On the other hand, we can not eliminate also 
an alternative interpretation which may be convention- 
ally denominated as the temporal one, i.e. genera- 
tion of some temporal wave processes in the course 
of the cosmological evolution (e.g 



Aref'eva & Koshelev 2008 Hirano & Komiya 



Morikawa 



1991 



2010). 



Although, it seems to be more difficult to explain a 
temporal structure formed by more-than-one periodi- 
cal processes. Note, however, that one temporal inter- 
val AT C = AD c /c = (330 ± 30) h" 1 Myr correspond- 
ing to the main quasi-periodicity indicated above could 
be consistent with the results of lAref'eva &: Koshelevl 
p008| . 

Let us note that results of our analysis could be quite 
sensitive to possible radial incompleteness of the DR7- 
Full sample. It is likely that some fainter LRGs at larger 
redshifts were missed by the survey. To minimize such a 



selection Kazin et al. ( 2010a ) restricted their analysis by 
redshifts z ~ 0.36 corresponding to the quasi-volume- 
limited sample. They performed additionally some an- 
gular and radial weighting procedures and showed that 
their results did not variate essentially. Using the same 



catalog Blake et al. (2011) extended similar analysis to 



higher redshifts z — 0.44 and obtained the close peak 
position in the correlation function. Here we analyse 
the whole DR7-Full sample of LRGs up to z max = 0.47. 
However, as it seen in Fig. [2] the distribution function 
N(??) is still quite representative at the edge of highest rj. 
All that may evidence in favour of an assumption that 
at least the scales ~ 100 hr 1 Mpc are not distorted 
strongly by the radial-selection effects. 

Fig. [7] gives examples of characteristic scales which 
have been obtained in the last decade by a few groups 
of authors as a result of analysis of large-scale distri- 
butions of galaxies and galaxy clusters. A part of data 
displayed in Fig. [7] is taken from Table 1 of |Ryabinkov| 
& Kaminker (20111). Additional segments of scales are 



introduced with using results of Hiitsi (2006) on the 



BAO incorporated in the power spectrum of the SDSS 



DR4 LRG sample, and Percival et al. (20101 on the 



disclosure of the BAO in the power spectra of both 
SDSS DR7 samples: LRGs and MGS, with inclusion of 
2dF Galaxy Redshift Survey data. We add also the re- 



cent results by Arnalte-Mur et al. (2011 1 on the wavelet 
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analysis of acoustic wave features in the spatial galaxy 
distribution. 

Fig. [7] should be considered only as an illustration 
and it includes rather nonuniform results. For in- 
stance, | Tago et aL ( 2002 ) yielded a characteristic pe- 
riod, (115±7) Mpc, of regular spatial oscillations of 
the two-point correlation function calculated for galaxy 
superclusters (see also Einasto et al. 1997a|b I. Two 
scale intervals chosen from Table 1 of Doroshkevich 



et al. (20041 correspond to the mean separations be- 
tween walls - the largest elements of the Large Scale 
Structure (LSS) in the galaxy distribution. The scale 
(74 ± 17) h" 1 Mpc was obtained on the basis of ra- 
dial (line-of-sight) measurements with using data of 
SDSS DR1 catalog (at z < 0.14), while the scale 
(60 ± 10) h^ 1 Mpc corresponds to the Las Campanas 
Redshift Survey data (at z < 0.11). Three significant 
periods, (127 ± 21) hr 1 Mpc, (73.4 ± 5.8 ) hr 1 Mpc, 
and (31. 7± 1.8) hr 1 Mpc, were revealed by Hartnett & 



Hirano ( 2008 ) in their Fourier analysis of radial galaxy 



distribution on the basis of SDSS DR5 and 2dF GRS 
data. The periodicity found by Ryabinkov & Ka minker 



(2011) in the radial distribution of QSO ALSs is also 



included. 

In spite of heterogeneity of the results, the scales 
in Fig. [7] may be formally subdivided into three 
groups: two scales belong to the interval (30 ± 
3) h" 1 Mpc, four scales belong to the wider interval 
(70 ± 21) h^ 1 Mpc, and the most representative inter- 
val is (121 ±26) hr 1 Mpc. Four periodical components 
discussed in the present paper get into these intervals; 
two most significant periods of our analysis belong to 
the third group. Especially good agreement of the main 



scale , (101 ± 6) h 1 Mpc, occurs with results obtained 



al. (2010), and Hiitsi 



by Blake et al. (2011), Kazin et al. (2010a), Percival et 



(2006) for the BAO scale. 



Actually, BAOs are discussed in literature as a scale 
of the sound horizon at the recombination epoch (e.g. 



Blake & Glazebrook 2003 Percival et al. 2007M dis 



playing itself as a single feature in the spatial correla- 
tion function (e.g. Eisenstein et al. 2005 Kazin et al. 



2010a Blake et al. 2011). This feature appears also 



as a series of regular variations (oscillations) imprinted 
in the power spectrum calculated for the spatial 3D- 
distribution of galaxies (e.g. |Eisenstein & Hu||1998 



et al 



2010 



Eisenstein, Hu fc Tegmark||1998| |Hiitsi||2006[ |Perciva 



Ross et al. ||2012[ ). They principally dif- 



fer from the quasi-regular variations of the radial (ID) 
distributions in real space discussed by Ryabinkov fc 
Kaminker (2011) and the present paper. However, it 



is not unlikely that both types of variations could be 
reconciled. 

The radial distribution variations probably bring out 
some spatial partly-ordered structure of matter in the 
early Universe (e.g. Einasto et al.||1997b| |Einasto et al. 



2011[ ) characterized by a long-range (or intermediate 
range) order. In particular, one can admit the possi 
bility that primordial acoustic perturbations, responsi- 
ble for the BAO, could carry traces of a partly order- 
ing formed at some early epochs (e.g. radiation-matter 
equipartition or recombination). In that case the ra- 
dial (ID) distributions, as well as the two-point 3D- 
space correlation function, would display a set of quasi- 
periodic features, which might reveal itself into a com- 
plex set of features (e.g. peaks) in appropriate power 
spectra. Note that possible presence of the second fea- 
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ture in the two-point correlation function at a scale ap- 
proximately double to the BAO scale ~ 200 hr 1 Mpc 



has been revealed recently by Ross et al. (2012); see 



also Martinez ct al. (2009) 



In any case, the existence of the large-scale periodici- 
ties, as well as the space-ordering hypothesis, needs fur- 
ther verifications based on statistical properties of dif- 
ferent cosmological objects over wider redshift regions. 
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A Effects of trend uncertainties 

Fig. [8] demonstrates variations of the power spectra 
obtained with using different procedures of the radial 
smoothed function (trend) subtraction (see Sect. [3]). 
We examine 105 831 LRGs from DR7-Full sample rep- 
resented in Table 1 by Kazin et al. ( 2010a[ ) . The upper 
panel in Fig. [8] is designed similar to Fig. [2] but with us- 
ing three different ways of the trend determination. In 
all cases the radial (line-of-sight) distribution function 
N(7j) is calculated for independent bins 6 V = 0.00033 
or c • 5 v /Ho — 1 hr 1 Mpc providing narrow spike- like 
variations of N(t^). The first trend function N tr (r?) is 
calculated by the least-squares method with using a set 
of parabolas (thick solid line) as a regression function 
for N(?7). The second one is calculated as a sum of the 
first five reciprocal Fourier harmonics, m=l, 5 (thick 
dashed line), with full filtration of all the rest harmon- 
ics 6 < m < 40. The third trend is calculated with 



using the radial selection function n(z) by Kazin et al. 
( |2010a[ ) presented on the World Wide Web (see Sect.[lJ) 



We rescale n(z) into N tr (77) and applicate a spline inter- 
polation for additional smoothing of the trend function 
(dotted line). Let us emphasize that x 2 -criterion does 
not admit to prefer one of three trends although the 
results of power spectrum calculations for each of them 
are qualitatively different. 

The lower panel displays the respective power spec- 
tra calculated with Eqs. Q and ^ and depicted by 
the same type of lines. One can see how the two most 
significant peaks at k = 0.062 h Mpc -1 (m=8) and 
k = 0.047 h Mpc" 1 (m=6) strongly depend on a trend 
obtained in one of the ways indicated above: the fea- 
tures may variate from strong double peak (dashed line) 
up to its cancellation (dotted line). Whereas for k > 0.1, 
i.e. for shorter waves, the effects of different trends 
are less essential. One can notice that in the upper 
panel the spline interpolation of the trend (dots) dis- 
plays long-wave oscillations and thus compensates the 
corresponding Fourier harmonics. To some degree such 
a compensation may occur for different trends subject 
to their sufficient complexity. The oscillations inherent 
in trends can have also as physical as nonphysical origin 
and in their turn need special statistical considerations. 

The indicated uncertainties make the binning ap- 
proach, including a procedure of the trend subtrac- 
tion, quite ambiguous at least in the cases of complex 
trends. Therefore, one needs a verification by a few 
statistical tests additional to the binning approach dis- 
cussed in Sect. [3j including the point-like treatment, 
to make more robust conclusions concerning reality of 
quasi-periodicities in the distribution of matter. 
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Fig. 8 Effects of a procedure of smoothed function (trend) 
determination on appropriate power spectra calculated in 
the binning approach. Upper panel: same radial distribution 
function N(rj) as in Fig. [5] but with three different trends: 
solid line - set of parabolas (same as N tr (??) in Fig. [5| , dashes 
- reciprocal Fourier transform employing first five harmon- 
ics, dots - radial selection function obtained from a comov- 
ing number density n(z) of Kazin et al. (2010a) with car- 
rying out a spline interpolation. Lower panel: three power 
spectra calculated with using Eqs. [4] [5] and the trends dis- 
played in the upper panel by the respective type of lines. 



